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Abstract: We study non- linear massive gravity in the spherically symmetric context. Our main 
motivation is to investigate the effect of helicity-0 mode which remains elusive after analysis of 
cosmological perturbation around an open Friedmann-Lemaitre-Robertson- Walker (FLRW) uni- 
verse. The non-linear form of the effective energy-momentum tensor stemming from the mass 
term is derived for the spherically symmetric case. Only in the special case where the area of the 
two sphere is not deviated away from the FLRW universe, the effective energy momentum tensor 
becomes completely the same as that of cosmological constant. This opens a window for discrim- 
inating the non-linear massive gravity from general relativity (GR). Indeed, by further solving 
these spherically symmetric gravitational equations of motion in vacuum to the linear order, we 
obtain a solution which has an arbitrary time-dependent parameter. In GR, this parameter is a 
constant and corresponds to the mass of a star. Our result means that Birkhoff's theorem no 
longer holds in the non-linear massive gravity and suggests that energy can probably be emitted 
superluminously (with infinite speed) on the self-accelerating background by the helicity-0 mode, 
which could be a potential plague of this theory. 
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1 Introduction 

The theory of massive gravity (see [1] for a recent review) has a very long history dating back to 
1939 when Fierz and Pauli [2] first wrote down the action describing a free massive spin-2 particle. 
Since then the progress in this area has been made sporadically and over the past few years this 
topic has been in a renaissance again because of the recent consistent construction of the non-linear 
massive gravity theory [3, 4]. 

The quadratic form of the massive spin-2 action is uniquely fixed as the linearized general 
relativity with the Fierz-Pauli mass term [5] by the ghost-free and tachyon-free conditions. Nev- 
ertheless, this linear massive spin-2 theory suffers from the van Dam-Veltman-Zakharov (vDVZ) 
discontinuity which gives, in the massless limit, the prediction for light bending of 25 percent differ- 
ence from that of General Relativity (GR) and hence is inconsistent with observational data. While 
a massive spin-2 particle has five physical degrees of freedom, namely the helicity-2, helicity-1 and 
helicity-0 modes, a massless spin-2 particle has only two helicity-2 modes. The vDVZ discontinuity 
stems from the fact that helicity-0 mode does not decouple with matters even in the massless limit. 

Non-linear effects are expected to relieve the inconsistency problem. Indeed, the linear analysis 
is no longer valid inside a certain distance from the source which is called Vainshtein radius. 
Because the Vainshtein radius goes to infinity as the graviton mass approaches to zero, we can 
no longer trust the above linear analysis in the massless limit. This opens the possibility that 
non-linearity may screen the helicity-0 mode at a scale rendering the theory compatible with 
observations. 

In compensation for screening the helicity-0 mode, the non-linearity gives rise to another 
obstacle, the Boulware-Deser ghost [6]. Non-linear interactions change the structure of the con- 
straints which in the linear theory eliminates unwanted degrees of freedom generally, hence ghost 
mode generally appears. In addition, the mass of the ghost is typically the same order as that of 
the massive graviton. We know from the experimental and observational data that the mass of 
graviton should be small, and thus we cannot simply ignore the ghost mode. The ghost issue has 
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been the core problem of massive gravity for a very long time and has been actively studied in the 
past few years [7-11] from the effective field theory approach [12]. 

Recently a two parameter non-linear massive gravity which is ghost-free has been developed 
by de Rham, Gabadadze and Tolley (dRGT) [3, 4]. In this construction, coefficients of higher 
order terms are chosen at each order such that in the decoupling limit the degree of freedom of 
the BD ghost is removed. Because of its construction, this theory is potentially free of BD ghost. 
Moreover, it was exactly proved that in this theory BD ghost does not show up at fully nonlinear 
order [13-16]. As a result, the form of the nonlinear terms are constrained and the mass and 
interaction terms are described with only three parameters. 

Once a gravitational theory is constructed, it is always interesting to survey its cosmological 
implication [17-31]. One of the phenomenologically interesting aspects is that small graviton mass 
naturally explains the smallness of the cosmological constant which we have observed. The basic 
idea is that due to the mass, a Yukawa suppression e~ mr comes in and the gravitational force is 
weakened at large scale. So even if we have a large underlying cosmological constant, say the one 
predicted by usual quantum field theory, it is effectively small in terms of gravity. The accelerating 
Friedmann-Lemaitre-Robertson- Walker (FLRW) solution of dRGT massive gravity was found in 
[17] where the spatial geometry can neither be flat nor closed [19] but open. 

The perturbation theory with general matter around self-accelerating universe solutions was 
studied in [18]. It was found that the kinetic terms of the helicity-0 and helicity-1 modes vanish, 
while the helicity-2 modes have time- varying mass. This is rather unexpected since usually massive 
gravity has five propagating degrees of freedom. This also makes it difficult to compare this non- 
linear massive gravity with GR because the structure of the propagations in the linearized theory 
of the former is basically the same as that in the linearized theory of the latter except for the 
helicity-2 part. 

The purposes of this paper are to study the difference between non-linear massive gravity and 
GR in open FLRW background, and to investigate the effect of helicity-0 mode. As a first step, 
we consider a spherically symmetric configuration in the vacuum, say outside of a star, but which 
asymptotically approaches to self-accelerating open FLRW background solution found in [17]. The 
corresponding solution in the framework of GR is the de Sitter-Schwarzschild solution, which has 
only one constant parameter, namely the mass of the star. We obtain the result that for this 
spherically symmetric setup, unlike in GR, an arbitrary time-dependent function appears in the 
solution in linear analysis. 

The rest of this paper is organized as follows. In Sec. 2 we review the non-linear massive gravity 
model and its self-accelerating solution. In Sec. 3 we analyze the spherical symmetric solutions 
whose asymptotic structure are the same as the self-accelerating solution. Sec. 4 is devoted to the 
summary of this paper and discussion. In Appendix A, we show the calculation related to the mass 
terms. In Appendix B, the derivation of the perturbative contribution from the Einstein-Hilbert 
term is shown. 

2 Review of dRGT massive gravity and its self-accelerating solution 

In order to be self-contained and to set up the notations, we briefly review some basic notions of 
massive gravity theory, the structure of the non-linear massive gravity proposed by de Rham et 
al. [3, 4] and the self-accelerating solution found in [17]. 

One important concept in massive gravity is that a priori background is needed in order to 
construct the theory. This seems to be a drawback as compared to GR which is background- 
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independent and can be described as geometry. However, one should not forget that by starting 
from the theory of free massless spin-2 particle propagating in Minkowski space and consistently 
adding non-trivial interactions one can arrive Einstein's GR. (For introduction to GR from this 
point of view, see chapter three of [32].) The background- independence of GR is a result of the 
consistent massless spin-2 theory. In this sense, the introduction of a priori background in the 
massive spin-2 theory is not intentional. 

From a practical point of view, since adding a mass term by only using the physical metric g^ 
is impossible, it is inevitable that another metric is required for massive gravity [1]. Particularly, 
the trace of g^ with itself only gives a constant. Hence, in order to introduce square or higher 
order terms, one requires another metric to do contractions and to take traces. The simplest way 
is to introduce a flat absolute background r\av Introducing spin-2 particle hnv on the flat absolute 
background, the physical metric g^ v is defined as 

and we assume that all matters feel the physical metric g^y. 

Note that hnv transforms non-covariantly under diffeomorphism and that the diffeomorphism 
symmetry is broken by introducing the non-dynamical absolute background. Nonetheless, we can 
construct an equivalent but covariant theory by introducing four scalar fields (fi a (a=0,l,2,3) called 
Stiickelberg fields. The relation between the Stiickelberg fields and the physical metric is defined 
as 



9yev — + H^ u , (2-2) 

where 

V = Vabd^ncb b (2.3) 

is called the fiducial metric. If we choose the gauge condition where <p a = x a , the fiducial metric 
becomes rj^ u and the field H^ u is reduced to hav- Note that since the Stiickelberg fields transform as 
scalars under diffeomorphism transformation, the fiducial metric Z^ v and the field H^ v transform 
covariantly. One replaces all occurrences of h^ u in the action with H^ u = g^ u — Z^ u , and then 
arrives in the final expression where all quantities transform covariantly and the diffeomorphism 
symmetry is restored. 

The action of the non-linear massive gravity proposed in [4] consists of two parts, the usual 
Einstein-Hilbert term Seh and a graviton mass term S mass : 



M 2 

Seh = ^t I dx A ^g-R, (2.5) 



S g = Seh + S mass , (2.4) 

K 

2 

S ma ss = M 2 Pl m 2 j d 4 xy/^g(C 2 + a 3 £ 3 + a 4 £ 4 ), (2.6) 

where R is the Ricci scalar constructed with physical metric g^ y . In preparation for constructing 
the mass term S maS s, we define a tensor fcf^v as 

K% = n - Wt, (2.7) 
= g^VabOa^d^ = Z% = cf"Z pv , (2.8) 
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where /i, u = 0, ■ ■ ■ ,3 and a, b = 0, • • • ,3. In order to avoid the appearance of BD ghost, the mass 
term S mass is constructed as [4] 

C2 = \([K] 2 -[K 2 ]), (2.9) 
/: 3 = ^([/C] 3 -3[/C][/C 2 ]+2[/C 3 ]), (2.10) 
£ 4 = 1 ([/C] 4 - 6 [/C 2 ] [/C] 2 + 3[/C 2 ] 2 + 8[/C] [/C 3 ] -6[/C 4 ]) , (2.11) 

where the square brackets denote trace operation, for example [IC] = ZC M ^. Indices are raised and 
lowered with respect to the physical metric g^, unless specified otherwise. In terms of W^ v and 
Z^u, £-2-, £3 and £4 are expressed as 

C 2 = Q - 2,[W\ + \[Wf - \[Z], (2.12) 

£3 = 4- 3[W] + [W] 2 - ±[W] 3 - [Z] + \[W][Z] - l -[WZ], (2.13) 

C 4 = l-[W] + \[W) 2 -\[Wf + ^[Wf 

-\[Z] + \[W][Z] - \[W) 2 [Z) + l -[Z) 2 - l -[WZ) + \[W][WZ] - \[Z\ (2.14) 

The gravitational equation can be obtained by the variation of the action with respect to the 
metric; 

-^G^iv + ^- — S mass = 0, (2.15) 

where G^ u is the Einstein tensor, id est, G^ u = R^ u — (R/2)g^ u . 

According to paper [17], we show the self-accelerating open FLRW solution in this theory. 
With the aim of finding FLRW solution, we fix the gauge such that the Stiickelberg fields are in 
the form of open chart of Minkowski spacetime : 



f(tWl + \K\(x* + y 2 + z 2 ), (2.16) 



■ AK\f(t)x, (2.17) 



'(0) 
;,2 



v (0 ) " V\K\f(t)y, (2.18) 

0?o) = VW\f(t)z. (2.19) 

In this gauge the fiducial metric respects the symmetry of the open FLRW space time: 

V = -iff^u + mfn^si, (2.20) 

where 

„ , , „■ ,9 ,9 ,9 \K\ (xdx + ydy + zdz) 2 

Qijdx'dx 1 = dx 2 + dy 2 + dz 2 - ' IV ,^,, / - 9 4-. 2.21 

1 + |ii |(ar + 2r 4- 2^) 

As for the physical metric, one considers an open {K < 0) FLRW ansatz; 

g^dx^dx" = -dt 2 + a 2 {t)Vl ij dx i dx : ' . (2.22) 
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Substituting (2.16)-(2.22) into the action and varying with respect to f(t) one obtains the solution 
for the Stiickelberg fields, which has two branches of solution, 



a _ 1 + 2a 3 + a 4 + ±\A + a 3 + aj - a 4 , . 

When a 3 and 04 are the orders of a small parameter e, we have 

X + = ^^+ to3 + C V ( e ) + ---> (2.24) 
a 3 + a 4 2(a 3 + a 4 ) 

X. = | + 0( e ) + ... . (2.25) 

The Friedmann equation obtained from this theory is basically that obtained by GR but with 
cosmological constant replaced by the graviton mass. In particular, 



(2.26) 



where 



2 A eff \K 

a = — H 

3 a 

£ = -^r, (2-27) 



A cfr = c±m2, (2.28) 



1 + a 3 ± a/1 + a 3 + a| - a 4 
(a 3 + a 4 ) 2 

In the small parameter limit (a 3 ,a 4 = 0(e) ), c± becomes 



c± = -± ■ ^ +a ")2 — " ~ ( 1 + "i - 2a 4 ± (1 + a 3 ) ^/l + a 3 + a| - a 4 ) . (2.29) 



c + = -t—^ ~ ' (Q3 = ^ - 3 , ( f ° 3 " ^ + O(0 + - , (2.30) 
(a 3 + a 4 ) 2 (a 3 + a 4 ) 2 4(a 3 + a 4 ) 2 

c_ = ^ + 0(e). (2.31) 

Note that the "positive" branch of solution becomes singular when a 3 and q 4 go to zero, while 
the "negative" branch remains regular. 

The Friedmann equations of nondinear massive gravity are the same as GR with cosmological 
constant, where the mass of the graviton plays the role of effective cosmological constant. This 
opens a window for explaining the smallness of the cosmological constant. Although we still need 
the small parameter, but the physical meaning is different. If we regard vacuum energy as the 
origin of cosmological constant, naive approximation using quantum field theory is of 10 120 order 
difference from observational results. Apparently, our knowledge about vacuum is very far from 
complete. On the other hand, because of the Yukawa suppression, rendering graviton small mass 
naturally explains the accelerated expansion no matter how large the vacuum energy density is. 
In fact, one cannot distinguish the contribution of the graviton mass and cosmological constant 
at the background evolution level. Specifically, if we add a cosmological constant term in the 
non-massive gravity action, one will arrive at Friedmann equations with a cosmological constant 
in addition to the graviton mass term which plays the same role as the former. Of course, this 
does not solve the whole problem but gives another physical interpretation. The quest for the true 
nature of vacuum still remains to be a very important problem in physics. 
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3 Spherically symmetric analysis around self-accelerating Background 

It was found in [18] that for open FLRW background the helicity-0 and helicity-1 modes have 
vanishing kinetic terms. Hence, in the linear analysis, the structure of the propagations of the 
cosmological perturbation theory on the accelerating background of dRGT massive gravity is the 
same as that of GR except that the helicity-2 mode in the former has a time-varying mass. Some 
natural questions then arise: Are there other observational consequences where we can see the 
contribution of helicity-0 mode? How about the helicity-1 mode? What will we find if we go 
beyond linear analysis to non-linear regime? In this work we will try to address the first and the 
third questions. 

Since in dRGT massive gravity five degrees of freedom can propagate on Minkowski back- 
ground, the situation that the helicity-0 and helicity-1 modes do not propagate on an open FLRW 
background at linear regime is very likely to change if we consider a background that is slightly 
different from open FLRW. One example is considering a star in an open FLRW universe. At a 
certain distance from the star, where linear perturbation is available but not far enough to ne- 
glect its gravitational effect, the helicity-0 mode propagates hence gives an additional force. The 
spacetime would no longer be de Sitter-Schwarzschild solution, as in the same configuration but 
with GR as the gravitational theory. Investigating this kind of spherically symmetric case which 
asymptotically approaches to open FLRW solution can not only give a strong phenomenological 
constraint on dRGT massive gravity, but also give us a better understanding on the helicity-0 
mode. 

A spherically symmetric metric can be generally written as 

ds 2 = _ e 2Z dt 2 + ^ w + 2 + a 2 r 2 e 2E dQ 2 ( g -Q 

1 — Kr K ' 

where f3 and E are the functions of t and r, 

dnf 2) =d9 2 + sin 2 9 dip 2 , (3.2) 

and for Q = '$? = f3 = E = this metric becomes exactly the open FLRW metric (2.22) which is 
obtained in [17]. Throughout the calculation we will fix the gauge such that the fiducial metric 
remain diagonal form to all order, 

Z^dxV = -(f(t)) 2 dt 2 + + \K\(f(t)) 2 r 2 dn 2 2) , (3.3) 

where / is defined in eq.(2.23). Gravitational equations of motion can be obtained by variations 
with respect to g^ v '. Although variations with respect to Stiickelberg fields give equations of 
motion, all of them can be reduced to the equations of motion from variations with respect to 
g^ u [33]. We show the contributions of the gravitational equations of motion from the mass term 
and the Einstein-Hilbert term separately. After that, we combine them and solve the linearized 
equations. 
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3.1 Contributions from the mass term 

In this subsection, we show the contributions from the mass term. The variations of [W], [WZ] 
and [Z 2 ] with respect to g^ v can be written as 



§p = W ^ = m^^ v + 2Wt^ + \ X±e ~ 9 ^ (3 ' 4) 

6[WZ] - (WZ) UV = + T^r K " V + lxle-* E (9»v - fa), (3-5) 



- 2Z, ia Z a t , = -2 det(Z)^ + 2 (W^ - 2Vditz) + 2X%e~^ E ( 9iiv - fa), (3.6) 



where 

det(Z) = (</V - g 0r g 0r )Z oZrr, (3.7) 



W T = \Jg 00 Z 00 + <T Zrr + 2Vdet(Z) , (3.8) 

are the determinant of the upper-left 2x2 part of Z^ v and the trace of the upper-left 2x2 matrix 
of W^y respectively, the upper-left 2x2 parts of g^ iu and Z^ are the same as those of g^ v and 
Z^y respectively, and the other components of g^ u and Z^ v are zero 1 . Namely, by using the matrix 
I a u = diag(l, 1,0, 0), they are exactly written as fa = g^ a I a v and Z^ = Z m l a v . 

Varying the £ 2 , £3 and £4 parts of the action with respect to g^ v yields the followings 



1 



1 5_ 



-gSg 



J d 4 x^g-£ 2 ^j = ~£ 2 g^ + (-3 + [W]) - X -Z yw , (3.9) 
d A xyf^g£i 



+ ("3 + 2[W] - l -[W} 2 + l -[Z]\ W, u + (-1 + \[W]\ Z, u - l - , 

(3.10) 



d x\J —g£^ 



\ug^ + (-1 + [W] - l -[Wf + l -[Wf + l -[Z] - \[W}[Z] + l -[WZ] ) 
1 1 1 „ 1, A _ / 1 1 \ 



-g $g 



+ ( -2 + 2^1" + - A \Z) J V + ( -3 + (WZ)^ " l Z ^Zl. (3.11) 

In order to investigate the difference between non-linear massive gravity and GR, we regard the 
mass terms in nonlinear massive gravity as effective energy-momentum sources in GR. From this 
point of view, (3.9), (3.10) and (3.11) are, then, related to the effective energy-momentum tensor 
as 

rpeS _ 2 5S mass . . 



1 Here, our coordinates are a; M = (t, r, 6, ip) 



-7- 



As calculated in Appendix A, each component of is written as 

-3 + 3A±e~ " 



T± = -2M 2 Pl m 2 q 



+a 3 (-2 + 3X ± e~ E - X±e~ 2E ) - y (l - X ± e~ E ) 2 } 



I V 2 



X±e~ 



+ a 3 



V-{l-X ± e- E y 



+ { (-3 + 2X±e- E ) + a 3 (-3 + 4X±e~ E - X 2 ± e~ 2E ) - q 4 (l - ^±e _jB ) 2 } W mn J ,(3.13) 
{-3 + 3 -X ± e~ E + a 3 (-2 + \x ± e~ E \ + f (-1 + X ± e~ E ) } 
+W r 1 1 - \x ± e~ E + a 3 Q - A ± e" s ) + f (l - A ± e" s ) } 



+ v /det(Z) 



-i + -x ±e - 



+ ^(-l + X ± e- E ) 



9ij, (3-14) 



where m,n = t,r and z, j = B,ip and the other components are zero. 

Interestingly, for E = 0, by using eq. (2.23) is reduced to -M Pl A cE9mn, namely there is 
no effective contribution except for the cosmological constant. Then, the form of can be fixed 
by the conservation law. The conservation law is obtained by take the divergence of eq.(2.15); 



■eff 



0. 



(3.15) 



While for v = i (i = 9,ip) the conservation law becomes trivial because of the symmetry, for v = t 

he relation between T^ n and Tf^ 
and v = r, it becomes the same constraint equation; 



and v = r it gives the relation between T^ n and T^. When T^ n = 0, in the both case for v = t 



0=(T l:j \ E= o + M 2 Pl A cS g l , j )5 ij 
= -2M 2 Pl m 2 [-3 + X± + a 3 (-3 + 2X±) + a 4 (-1 + X±)] [ X± - W T + 



v /det(Z) 



(3.16) 



where we use eq.(2.23) and (2.28). Since is proportional to gij, this obviously results in 

Tij\ E =o = -M Pl A cmj . (3.17) 

Therefore, for E = 0, the contribution from the mass term is completely equal to the cosmological 
constant and its solution must be the same as the gauge fixed one of GR with the cosmological 
constant, i.e., de Sitter-Schwarzschild with the following gauge choice 2 : 



E = 0, 



if 



l-Kr 7 



■2<P 



a 

W\ 



-2<I< 



1) 



(3.18) 
(3.19) 



We linearly expand T~ since in subsection 3.3 we analyze the linearized equations. Assuming 



$, /3, ^ and E are all the same order as a small parameter e and linearly expanding (3.9), Tj£ can 



With the different gauge choice, the de Sitter-Schwarzschild solution was obtained in [20-23]. 
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be expanded as 



T r f + M' Pl A eS g r 



-2M 2 Pl m 2 g C ± 1 



\K\ 



Eg rr , 



(3.20) 



T$ + M 2 Pl A eS9ij 



-M 2 Pl m 2 g C ± 1 



(3.21) 
(3.22) 



C± = {(-3 + X±) + a 3 (-3 + 2X±) + a 4 (-1 + X±)} X±, 
and the other components of + M pl A e sg ^J} are zero. 
3.2 Contribution from the Einstein-Hilbert term 

In this subsection, we show the linear contribution from the Einstein-Hilbert term. Since the form 
of the Einstein-Hilbert term is the same as that in GR, the linearized parts of the contribution 
from the Einstein-Hilbert term must be written with the gauge invariant combinations of GR. 
Therefore, it is convenient to define new variables A and B as 



A = ^>- 



a 2 r 2 
1 - Kr 2 
H 
H 2 



H 1 $ 



1 - Kr 2 
Kr 2 \ 
1 - Kr 2 J 



rH/3 + 



a 2 r 2 



H 



:E' 



Kr 2 ) 
E 



HE 



E 



1 — Kr 2 



rE', (3.23) 



H(l-Kr 2 ) H 2 



H . H E 
+ --^rE' + 



H 2 1 - Kr 2 ' 

(3.24) 



where prime and dot denote derivative with respect to radial and time coordinate respectively, and 
which are gauge invariant variables in GR. (See Appendix B.) We stress that in massive gravity, 
since there are no gauge degrees of freedom, these variables A and B are not gauge invariant 
variables but are just the linear combinations of the physical variables /3, \E r and E. These 
variables are useful to write the linear contributions from the Einstein-Hilbert term and to discuss 
the difference between massive gravity and GR. 
Equipped with these two variables, we have 



-Rgoa + A c ff5oo : 
1 



M^^) A+ 2 {1 -K^ A , + 2H2B + 0(e2) 



a 2 r 



2 H 



-Rg 0r + A c s9or = - \-A + HB\+ O(e'), 



H' 



-Rg rr + A c gg r 



H 
H 2 



l\A 



H 



A-HB 



+ 0(e 2 ), 



(3.25) 



r(l - Kr 2 )B' - a 2 r 2 {HB + HB) - (2Kr + 3a 2 r 2 H 2 )B 



r(l-Kr 2 ) .. Kr 2 ■ , 9n 
'-A' + —^A - r(l - Kr 2 ) 



H 



H 



A' 



H 



+Kr'\l- w \A + (D{e 



a 2 r 2 ^ 1 ^ 



3.3 Linear analysis 

In order to investigate the properties of massive gravity theory, we focus on the week gravity 
case and solve the linearized equations. Instead of the four variables /?, and E, we use the 
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variables A, B, ^ and E as independent variables. It is useful because equations become simpler 
and because A and B are gauge invariant combinations in GR. The contributions of the linear 
order is obtained in § 3.1 and § 3.2 and the linearized equations of motion become 



-A + 



2(1 - Kr" 1 

a 2 r 2 



2(1 - Kr 

a 2 r 



0. 



-A' + 2H 2 B = 0, 



H 

IP 



l\A 



H 



A-HB 



-2m 2 a C± 



1-Kr 2 



aH 



E. 



r(l - Kr 2 )B' - a 2 r 2 (HB + HB) - (2Kr 2 + 3a 2 r 2 H 2 )B - r ^ Kr ^ A' 



H 



Kr 2 f 
+ ^ r A-r(l-Kr 2 )ll 



H 

H 2 



A' 



Recall that H = K/a 2 , hence (3.27) gives 

B 

Substituting this into (3.26) one obtains 

A' = 



K 



a 2 H 2 
(1 - 3Kr 2 



-m 2 C±a 2 r 2 1 



.4. 



aH 



■A. 



r{l-Kr 2 )' 

With these two relations, the last two equations of motion (3.28) and (3.29) become 



Jj2 r 2 

(1 - Kr 2 ) 
H 2 



H 2 A + HA 
H 2 A + HA 



-2m 2 C ± 



1 - Kr 2 



1 



-m z a C±a 2 r 2 I 1 



aH 



aH 
\f\K\, 

{^ + E) 



E, 



which yield 



* = -2E. 

Solving (3.31) and substituting the solution into eqs.(3.30), (3.32) and (3.34), we have 

1 M(t) 



A 
B 

E 



SirM 2 pi ar(l - Kr 2 ) ' 

1 K M(t) 
~8ttM 2 a 2 H 2 ar{l-Kr 2 )' 



1 



aH 



%irM 2 pi m 2 g C± \ 

1 2 

~MM 2 pi m 2 C± 



1 



K\ J 
aH 

fW\ 



M(t) 
a°r 

a 3 r 3 



(3.26) 
(3.27) 

(3.28) 



(V + E). (3.29) 



(3.30) 



(3.31) 

(3.32) 
(3.33) 



(3.34) 
(3.35) 

(3.36) 
(3.37) 



where M{t) is an arbitrary time-dependent function. The solution has one time-dependent func- 
tional parameter, which is very different from that in GR. Only when M{t) is a constant, the 
solution is the same as the linearized de Sitter-Schwarzschild solution and M{t) is corresponding 
to the mass of a star (or a black hole). Then, E and ^ are zero, which is consistent with the 
non- linear result in § (3.1). 
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4 Summary and Discussion 



In this paper, we have analyzed the spherically symmetric solutions with asymptotically open 
FLRW solution in non- linear massive gravity. In § 3.1, we have presented the non- linear form of 
the effective energy-momentum tensor which comes from the mass term. Only in the special case 
where E = 0, the effective energy-momentum tensor is completely the same as that of cosmological 
constant. Otherwise, other effects show up. This gives a possibility to discriminate the massive 
gravity theory from GR. 

In § 3.3, we have solved the gravitational equations of motion in the vacuum in linear order 
of metric perturbation and obtained the result that the solution depends on a time-dependent 
parameter. In GR, the parameter is a constant and corresponding to the mass of a star. While in 
the analysis of gravitational propagation [18] the helicity-0 and helicity-1 modes do not show up 
and the difference from GR in linear order appears only in the mass of helicity-2 modes, in our 
analysis we can see the large difference from GR even in linear order. This is because we analyze 
not the gravitational propagation but the gravitational potential which is singular at the center. 
Because of the central singularity, our result is not included in the analysis of the propagation. 

Our result means that in massive gravity Birkhoff's theorem does not hold true and sug- 
gests that energy can be gravitationally emitted even in spherically-symmetric configuration. The 
existence of the helicity-0 mode is expected from this result because the spherically-symmetric 
gravitational wave is probably related to the helicity-0 mode. This gives a hint for the question 
if the propagation of the helicity-0 mode accidentally disappears only on open FLRW solution. 
The fact that the parameter only depends on time but not on the radial coordinate implies su- 
perluminous behavior. In a construction of a concrete solution, we choose the parameter at some 
point on each time-constant hypersurface, such as the surface of a star, and then, metric at other 
points on the same time-constant hypersurface is fixed. This means that its information propa- 
gates to all points on the same time-constant hypersurface, namely the speed of the propagation 
is infinite. Infinite speed of propagation can be realized when the coefficient of the kinetic term 
(time-derivative term) in the action vanishes. Therefore, we can conclude that our result shows 
the sign of superluminous (infinite-speed) behavior of the helicity-0 mode in the self-accelerating 
background. 

If the matter is coupled with the superluminous helicity-0 mode, the infinite-speed propagation 
makes the situation hazardous, although we do not show the coupling of gravity with matter in 
this paper. The superluminous helicity-0 mode conveys the information of matter to infinity for 
an instant. It can be observed as the violation of the energy conservation. In order to avoid the 
pathology in massive gravity, decoupling between matters and the helicity-0 mode (and probably 
also helicity-1 modes) is needed. 

In order to see the disappearance of the kinetic term directly, it is better to analyze the 
perturbation on the de Sitter-Schwarzschild solution and to take the limit that the mass parameter 
of the solution goes to zero. If our expectation that the kinetic term of the helicity-0 mode 
accidentally disappears only on open FLRW solution is correct, the helicity-0 mode can propagate 
on the de Sitter-Schwarzschild background. Then if we take the massless limit of black hole mass, 
we can check the disappearance of the kinetic term. Perturbation on the de Sitter-Schwarzschild 
background will be discussed in an upcoming work. 
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A Calculation of contribution from mass term 

In this appendix, we show the brief calculations about the mass term. On the metric (3.1), the 
concrete forms of [W], [Z], [WZ] and [Z 2 ] become 

[W] = W T + X±2e~ E , (A.l) 

[Z] = W T - 2^det(Z) + X 2 ± 2e~ 2E , (A.2) 

[WZ] =W T - 3W r ^det(Z) + 2X\e~ ZE ', (A.3) 

(A.4) 



[Z 2 ] =W T - 4W TV / det(Z) + 2det(Z) + 2X^e~ iE \ 
Substituting them into eqs.(2.9)-(2.11), we have 

l -C 2 = 3 - 3X ± e~ E + \x\e- 2E - W T Q - X ± e~ E ^ + ^x/deT(Z), (A.5) 
\C 3 = 2- 3X ± e~ E + X 2 ± e~ 2E + [~\ + 2X ± e~ E - ^e" 2 ^ W T + (l - X ± e~ E ) V / deI(Z), 



\U = \{1-X ± e- E ) 2 



1 - W T + v /det(Z) 



(A.6) 
(A.7) 



Then, the variation of £2, £3 and £4 of the action with respect to g^ v can be written as 



I 5_ 

1 



d xsj—g£; 



£ 2 g^ + (-3 + 2X±e~ E + W T ) W y 



-3 + 3X±e 



-Xle 



2 -IE 



+ W T 



X±e 



9nu + (-3 + 2X ± e~ E ) W l 



[lis 



-3 + \x ± e~ E + 



3 

2 2 



\x±e-A Wt-^VMZ)) (9 



\c*9v, + ("3 + 4X ± e" s - Xle~ 2E + (2 - 2X±e- £ )W T - Vdet(Z)) W, 



+ ( -1 + X±e- + -W T J V - 3 OT P 



-2 + SXie"- 5 - x£e _2i5 + 



- - 2X ± e~* + -X^e 



2 -2E 



W T 



+ (-3 + AX±e~ E - X 2 ± e~ 2E ) + 



-2 + -X ± e~* + l--X ± e 



W T 



-1 + ^X±e- E ) Vdet(Z) 



(fiV^ fl'/if) 1 



(A.9) 



1 6 



= -\U9»V + f-1 + 2X ± e~ E - X 2 ± e~ 2E +W T (l - 2X±e~ E + X 2 ± e~ 2E ) 



+ (-1 + 2X±e~ E ) Vdet(Z)} 
+ j-i + X ± e~ E - l -X 2 ± e~ 2E + W T Q - X ± e- E ^j - V 

> - *±e- £ ) 2 ( ~ (1 - W T ) ^ - WvJ) 
+ (1 - X ±e - B ) (-1 + W T - Vdet(Z)) (e^ - ^) , (A.10) 



where g^ is defined under eq.(3.6) and VV M1/ = W^ a I". The effective energy-momentum tensor 
Tf^ can be written with the linear combination of eqs.(A.8)-(A.10) as eqs.(3.13) and (3.14) and 
linearly expanded as 

T c ± + M 2 Pl A cS g mn ~ -2M|,m^£: (x ±5mn , - WPg mn + ) , (A.ll) 
Tff + M| z A efT ^ ~ -M 2 Pl m 2 g C ± [ 1 - -i== J (E + *) ^, (A.12) 



where C± is defined in (3.22) Background values of Wt and W mn are, respectively, 



7(0) 




a 



+ 1 I , (A.13) 



From these, we know the only rr-component of linearized (T mn — Ag mn ) becomes non-zero: 

T r f + M Pl A eS g rr ~ -2M|,m2C± [ 1 - -^== J ^rr- (A.15) 



B Derivation of contribution from the Einstein-Hilbert term 

In this subsection, we derive the linear contribution from the Einstein-Hilbert term. Because the 
contribution from the Einstein-Hilbert term should be the same as that of GR, we can use the 
same technique as in GR. In concrete terms, we construct the gauge invariant variables and write 
the perturbative terms with them. Since the background we consider is accelerated by the effective 
energy momentum tensor, we add the same contribution in GR by hand, namely we consider the 
perturbation of 

Rnv - 2 R 9/j.v + A cS9iiu- (B.l) 

In order to keep spherical symmetry, we would only consider diffeomorphisms in the temporal and 
radial coordinate: 

x» -> x" 1 = x» + e, e = (£', f , 0, 0) = (C, £, 0, 0). (B.2) 
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Under such gauge transformation, the metric perturbations transform as 



Kr 



(B.3) 



E E - H( 

r 

where prime and dot denote derivative with respect to the radial and time coordinate respectively. 
Then, the gauge invariant variables A and B are written as eqs.(3.23) and (3.24). Equipped with 
these two variables, the perturbation of (B.l) can be expressed as eq.(3.25). 
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